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Executive Summary 



Water depths in the near-shore zone of western Isla de Vieques arc calculated 
using Landsat 4 Thematic Mappet imagery and two different algorithms. The 
first is a new algorithm that docs not require a sea-bottom-type clustering; 
the second is the dual-band-ratio algorithm. TM bands one and two (0.45- 
0.52 pm and 0.52-0.60 pm, respectively) arc water penetrating and provide 
capability for calculating water depths to about 20 m in clear water A 
comparison of results of remotely sensed bathymetry using J/* two different 
models is made. 
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A Comparison of Models for Remotely Sensed Bathymetry 


Introduction 

Remotely sensed bathymetry from multispectral 
imagery has been studied for several years. 1 ’ Three 
models have been devised to enable calculation of 
water depth, the so called single-band algorithm,* the 
two-band ratio method.' and the multiband method. 6 
The first two models have the disadvantage of requiring 
a bottom-type clustering to minimize the bottom 
reflectance variation. The third model requires many 
(>2> water-penetrating bands, because present sensors 
are predominandy optimtzed for land remote sensing 
and do not have sufficient water-penetrating bands to 
satisfy this model. 

The single-band algorithm and the dual-band ratio 
method have been tested with good results.’ 6 However, 
a new algorithm we call the linear mulliband method 
has recently been introduced 10 that does not require a 
bottom-type clustering. This new algorithm holds the 
promise of yielding absolute bathymetry in cases where 
the number of different bottom types is equal to or less 
than the number of water-penetrating bands in the 
imagery. 

We recently tested the linear multiband method using 
bands one and two of Landsat Thematic Mapper 
imagery in the vicinity of Isla de Vieques. The results 
are compared to bathymetric measurements made 
using the dual-band ratio method. The calibration data 
used in both methods were taken from a 1:25000 
combat chart of the region. With the linear multiband 
algorithm we achieved improved results without 
performing the time-consuming process of bottom- 
type clustering and classification required by the dual¬ 
band ratio model 

Background 

Remotely sensed bathymetry from multispectral 
imagery generally is based upon a simple reflectance 
model;* the radiance in wavelength band i at water 
depth r is given by 

L, =i„ + c t R ai ixpl-ik'Z), 


where Z., is the radiance value in band i, L tm is the aver¬ 
age signal over deep water, c, is a constant that is a 
function of several optical parameters, R at is the 
bottom reflectance in band t over bottom type a, and 
i, is the diffuse attenuation coefficient. 

Solving for z, one obtains the formula 

Z - LN (c,R m ) !2k l -X l l2k l . 

where we have adopted the convention that X , = 
LS(L r L „) 

This single-band reflectance model assumes that the 
bottom retlectance is constant over the bottom type, 
that the atmosphere and the sea state are uniform, and 
that other background optical effects are cither uniform 
or constant throughout the image. 

To reduce errors due to the variation of the bottom 
reflectances, a two-band ratio method (or dual-band 
method) was devised.’ In this algorithm, the depth is 
given by 

Z*(l’2<k ;-*;)) 'LS <c,Rai'Ct*ai) * X, - X,l ■ 

where the subscripts / and 2 indicate different bands. 
In this method, we assume that changes in the 
bottom reflectances occur in such a way that the ratio 
c iRqiIc^R ai remains constant. 

Another method, which we call the linear mulliband 
method, gives the depth by 

Z--la i a'2k i )aHlc l R a j -X.) 

where the sum is taken over several bands and the 
weights co t satisfy the constraint that Iw,= / Paredes 
and Spero, ’ generalizing the assumption that the ratio 
remains constant, assume there are con¬ 
stants i, and a independent of bottom type such that 

<c l R 0 i)' ; ’lc 2 Ra 2 ) i ‘(c l Ra>^ = 

(ciR^'IczUbzi^lc^bj)^ ■■■ = 

= a 
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With this assumption, they show that certain weights, 
co i can be found such that, when used in the above 
multiband linear depth equation, produce an equation 
for Z independent of the bottom reflectance and depend 
only upon the values s,: 

Z = Ul2£ c ki)(l-S l X l -S t X 1 -.}. 

Alternatively, this equation may be wriaen 

Z = Ac + A 2 X 2 + • •• + A t X f , 

where the coefficients A 6 .Ai . A , are constants inde¬ 

pendent of the bottom type at which the depth is being 
calculated. (See ref. 7 formone details.) 

It has been shown that the assumption underlying the 
two-band ratio method is not satisfied by Landsat TM 
data 1 and it has been suggested’ that tne linear multi¬ 
band algorithm, with the generalized ratio assumption, 
may be more accurate given more bands than bottom 
classes, as the coefficients in the algorithm do not 
depend on the bottom reflectance. In this paper we use 
the linear multiband algorithm on two channels of a 
Landsat 4 TM scene and demonstrate its improved per¬ 
formance over the two-band ratio method. 

Procedure 

The Landsat scene used in this study (scene number 
5032614162, taken January 21, 1985) contains Isla de 
Vieques, an island off the southeast coast of Puerto 
Rico in the Caitibean Sea. The water clarity in this area 
is extremely good, and bottom reflectances in the blue 
and green channels (TM bands one and two) are detect¬ 
able for depths up to 25 m. Bottom reflectances in the 
red channel (TM band three) are detectable for depths 
up to about 6 m; sensor noise (striping) over the water 
renders this band useless at greater depths. 

The general procedure for this comparison of the 
three methods was first to convert a National Aeronau¬ 
tics and Space Administration/Goddard Space Flight 
Center computer compatible tapes to an image data file 
and to georeference this data file against a Defense 
Mapping Agency 1:25000 combat chart. See reference 
8 for details on the georeferencing procedure. This 
combat chart contains not only bathymetric soundings, 
but also land feature identifications that aid in the geo- 
referencing process. A set of approximately 600 cali¬ 
bration points was then selected by recording depths 
from the chart; about 300 of the points were used in our 
regression fit and the remaining 300 were used as a test 
set to check the calculated depth against the actual 
depth. 

The average deep-water signal for each water- 
penetrating band (in this case channels 1 and 2, roughly 
corresponding to the blue-green and green portions of 


the visibly. spectrum, respectively) was calculated to 
obtain the L im values. 

Polygons were drawn on the image to outline the 
deep water (no visible ocean bottom) region in each 
band. The average intensity of the pixels contained in 
the polygon was used as the f.,_ values. These values 
are subtracted from the corresponding L value to adjust 
the signal values for atmospheric scattering, etc. 

A linear regression of the model equation against the 
(appropriate) calibration points is run, which produces 
coefficients and, hence, equations for the depth Z at any 
pixel. The test set of calibration points is used to test 
the fit of this regression. These equations are then used 
to produce a bathymetric image that can be processed 
further (smoothed, contoured, pseudocolored, etc.) 
depending upon the desired application. 

Results and Analysis 

The linear regression was run on the equation 

Z = Ag + A;X/ + A 2 X 2 

corresponding to a two-channel multiband linear 
model. For this image we determined L im = 68 and 
f. 2 _ = 17. The regression yielded the following values 
for the constants; A 0 *ll 2.A,*3 H.A 2 =-676 The 
multiple correlation coefficient was 0.86. 

Scatter plots indicated an apparent correlation of 
error with depth, so we plotted the test residuals of 
measured depth versus calibration depth for three dif¬ 
ferent depth ranges, 0-5 m. 5-10 m, and 10-16 m. 
Residual means and RMS values are given in Table 1 
It was seen that we consistently underestimated the 
depth in the shallower water and overestimated the 
depth in the deeper water. 

Using the test set of calibration point data, the model 
yielded an overall residual mean less than 0.2 m and an 
overall RMS less than 1.9 m. Performing the regression 
in the depth range of 0-5 m, the RMS error was 0.88 m; 
in the range 5 -10 m, the RMS error was 1.25 m; and in 
the range 10-16m, the RMS error was 1.00m. 

For comparison purposes, the dual-band ratio 
method was applied as well. To minimize error, the 
image was clustered to locate areas of similar bottom 
reflectance using a supervised statistical clustering 
routine and the maximum likelihood classifier. The 
algorithm was then regressed against calibration points 
in each cluster separately. An overall residual mean of 
0.093 m and an RMS error of 2.66 m was obtained 
(Table 2). In the depth range of 0 - 5 m, the RMS error 
was 0.94 m; in the range 5 - 10 m. the RMS error was 
1.56 m; and in the range 10- 15 m, the RMS error was 
1.86m. 
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Table t. Errors in depth calculation from 
linear multiband model. 


Depth Range 

Residual Mean 

RMS 

0 - 5 m 

0 002 m 

0.88 m 

5 - 10 m 

- 0.39 m 

1.25 m 

10 -15 m 

- 0 21 m 

1.00 m 

0 - 15 m 

0.20 m 

1 86 m 

Table 2. Errors in depth calculation from 

dual-band ratio model. 


Depth Range 

Residual Mean 

RMS 

0 - 5 m 

0.12 m 

0.94 m 

5 - 10 m 

- 0 15 m 

1.56 m 

10 -15 m 

- 0 25 m 

1.86 m 

0- 15 m 

0 10 m 

2 66 m 

In comparing the 

two methods, both algorithms 

underestimated the 

depth in shallow 

water and 

overestimated it in 

deeper water. This tendency 


produces a larger overall error than is obtained when 
considering the indicated depth ranges. It can be seen 
that the linear multiband method yields somewhat 
improved results, even when only two bands arc 
available. Moreover, this method does not require the 
clustering and classification routines to discern areas 
of similar bottom reflectance, and a coasidcrablc 
savings in CPU processing time (one to several hours 
on a VAX 11/780 per 512 X 512 pixel image) is 
realized. 

Applications 

A major advantage of calculating water depths from 
remotely sensed data is the reduction in preparation 
lime and in the cost of producing bathymetric chans. 
Moreover, the digital nature of the product lends itself 
to easy and rapid editing and updating. 

The extrapolation of calculated bathymetric 
coefficients from one ponion of a scene to another 
portion is easily accomplished due to the digital nature 
of the data. The authors intend to study the validity of 
such extrapolations in future work as this has the 
potential of being a major mapping, charting, and 
geodesy application. 
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